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Abstract 

Genetic switch systems with mutual repression of two transcription factors are studied using 
deterministic and stochastic methods. Numerous studies have concluded that cooperative binding 
is a necessary condition for the emergence of bistability in these systems. Here we show that for a 
range of biologically relevant conditions, a suitable combination of network structure and stochastic 
effects gives rise to bistability even without cooperative binding. 
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Recent advances in quantitative measurements of gene expression at the single-cell level 
have brought new insight on the importance of stochastic fluctuations in genetic cir- 



cuits j3[. Populations of genetically identical cells show variability due to fluctuations. The 
role of fluctuations is enhanced due to the discrete nature of the transcription factors and 
their binding sites, which may appear in low copy numbers 0, . Stochastic behavior may 
invoke oscillations 0, and spatio-temporal patterns 0, Q| , which are unaccounted for 
by macroscopic chemical rate equations. Genetic circuits with feedback mechanisms often 
exhibit bistability, namely, two distinct stable states which can be switched either spon- 
taneously or by an external signal U, 12, 3|- To qualify as a switch, the spontaneous 



switching rate must be much lower than the rates of the relevant processes in the cell, 
namely transcription, translation, binding and unbinding of transcription factors. In partic- 
ular, genetic switches such as the phage A switch, enable cells to adopt different fates jl^ . 
The toggle switch is a simple genetic circuit that consists of two proteins, A and B, with 
concentrations [A] and [B], respectively, which negatively regulate each other's synthesis (by 
concentration we mean the average copy number of proteins per cell). The production of 
protein A is negatively regulated by protein B, through binding of n copies of B to the A 
promoter (and vice versa). This process can be modeled by a Hill function, which reduces 
the production rate of A by a factor of 1 + A; [5]", where k is a parameter and n is the Hill 
coefficient In case that n = 1 the binding of a single protein is sufficient in order to 
perform the negative regulation, while for n > 1 the cooperative binding of two or more 
proteins is required. In numerous studies of the toggle switch system it was concluded that 
cooperative binding is a necessary condition for the emergence of the two distinct stable 

Q) 12, ll^ 13 • ^Iso observed that in presence of 



states characteristic of a switch 



cooperative binding, stochastic effects contribute to the broadening of the parameter range 
in which bistability appears |2o| . 

In this letter we show that stochastic effects enable bistability even without cooperative 
binding of the transcription factors to the operator, namely for Hill coefficient n = 1. 
Furthermore, bistability takes place even when the active proteins appear in high copy 
numbers. These results emphasize the necessity of stochastic methods in the analysis of 
genetic networks, even under conditions of high concentrations. 

The mutual repression circuit, referred to as the general switch is described by the 
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rate equations 

[A] = gA{l - [rB]) - dA[A] - ao[A] (1 - M) + 

[B] = gsil - M) - dB[B] - ao[B] (1 - [r^]) + ai[rs] 
[r'A] = ao[A]{l- [r-A]) - aiir^] 

[rB] = ao[B]{l-[rB])-ai[rB], (1) 

where gx (s~^), X = A, B, is the maximal production rate of protein X and dx (s~^) is 
its degradation rate. For simphcity, we ignore the mRNA level and take the processes of 
transcription and translation as a single step of synthesis [2^. The bound repressors are 
considered as separate species rx and their concentrations are given by [rx], providing much 
insight into the repression process [2^. Here, is a bound A protein that monitors the 
production of B, while is a bound B protein that monitors the production of A. Since 
there is a single promoter of each type, < [rx] < 1- The parameter (s~^) is the binding 
rate of proteins to the promoter and ai (s^^) is the dissociation rate. 

It is commonly assumed that the binding-unbinding processes are much faster than other 
processes in the circuit, namely ao, «i ^ dx,gx- This means that the relaxation times of 
[rx] are much shorter than other relaxation times in the circuit. Under this assumption, one 
can take the time derivatives of [rx] to zero, even if the system is away from steady state. 
This brings the rate equations to the standard Michaelis-Menten form 

[A] = gA/{l + k[B]) - dA[A] 

[B] = gB/{l + k[A])-dB[B], (2) 

where k = ao/ai is the repression strength. For a given population of free X repressors, the 
parameter k controls the value of [rx]- The limit of weak repression, [rx] <^ 1, is obtained 
when k[X] <^ 1, while the limit of strong repression, [rx] — 1, is obtained for k[X] ^ 1. 
These equations turn out to have one positive steady-state solution, thus at the level of 
rate equations this system does not exhibit bistability. For symmetric parameters, where 
Qa = Qb = g and dA = dB = d, this solution is [A] = [B] = [(1 + Akg/dY/"^ — 1 1 2k. 

In order to account for stochastic effects, the master equation approach 0, 2^ |2^ is ap- 
plied. In the master equation, the dynamic variables are the probabilities P{Na, Nb, rA, rB) 
for a cell to include Nx copies of free protein X and rx copies of the bound X repressor. 
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where Nx = 0, 1, 2, . . ., and rx = 0, 1. The master equation for the mutual repression circuit 
takes the form 



P{Na, Nb, Ta, Tb) = 9ASrs,oP{NA " 1, A^^, TA, Tb) + gBSr^,oP{NA, Nb - 1, Ta, Tb) 

+dA{NA + 1)P{Na + 1, Nb, ta, tb) + ^^(iVB + l)P(iV^, Nb + 1, r^, tb) 
-{9ASrs,o + 9BSrA,o)P{NA, Nb, va, Tb) - {cIaNa + dBNB)P{NA, Nb, va, vb) 
+ao[{NA + l)5r^,iP{NA + 1, Nb, 0, r^) + (Nb + l)5.^,iP(iVA, Nb + 1, r^, 0)] 
+ai[6r^,oP{NA-l,NB,l,rB) + Sr^,oP{NA,NB - l,rA,l)] 

-ao{NA6r^,o + NB6rs,o)P{NA, Nb, ta, tb) - + 6rs,i)P{NA, Nb, ta, tb), (3) 

where 6ij = 1 for z = j and otherwise. The gx terms account for the production of 
proteins. The dx terms account for the degradation of free proteins, while the ao (c^i) terms 
describe the binding (unbinding) of proteins to (from) the promoter site. The average copy 
numbers {X), where X = Na, NB,rA,rB, are given by {X) = J2 NP{Na, NB,rA,rB) where 
the sum is over all integer values of A^^ and A''^ up to a suitable cutoffs and over r^, = 0, 1. 
Note that for distributions that are skewed or exhibit several peaks, the average does not 
reflect the actual behavior in a single cell. 

To analyze the role of fluctuations in this circuit we have calculated the probability 
distribution P{Na,Nb) = '}2rA,rB ^^-^^^ ■^B,fA,fB)- We used the symmetric parameters 
g = 0.05 (s~^), which correspond to average production time of 20 seconds, and d = 0.005 
s~^) which means degradation time of 200 seconds, in agreement with experimental results 
2J| . To examine a broad range of relevant values of k we performed two sets of simulations. 
In the flrst set we chose ai = 0.5 (s~^) and varied ao, while in the second set we chose 
aQ = 0.5 (s^^) and varied ai. We conflrmed that the population of free proteins depends 
only on the ratio, k. 

Under conditions in which the promoter sites are empty most of the time, namely rx -C 1, 
the repression is weak and the steady state solution exhibits coexistence of A and B proteins 
in the cell. In this case the distribution P{Na, Nb) exhibits a single peak [Fig. ^a)]. In this 
case, the values of (Na) and (Nb) obtained from the master equation coincide with [A] and 
[B], obtained from the rate equations. For strong repression, the distribution P{Na, Nb) 
exhibits a peak in which the A population is suppressed and a peak in which the B population 
is suppressed, as expected for a bistable system. However, a third peak appears near the 
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origin, in which both populations of free proteins diminish [Fig. ^b)]. This peak represents 
a dead-lock situation, caused by the fact that both A and B repressors can be bound 
simultaneously, each bringing to a halt the production of the other specie. This result is in 
contrast to the rate equations which exhibit a single solution, [A] = [B] , for the entire range 
of parameters. Below, we present three biologically sensible variants of the circuit in which 
the third peak is suppressed, giving rise to a bistable switch. 

Consider the exclusive switch, where there is an overlap between the promoters of A and B 
and thus no room for both to be occupied simultaneously. Such a situation is encountered 
in nature, for example, in the lysis-lysogeny switch of phage A It was shown that 

in presence of cooperative binding, the exclusive switch is more stable than the general 
switch This is because in the exclusive switch the access of the minority specie to the 
promoter site is blocked by the dominant specie. Here we show that in the exclusive switch, 
stochastic effects give rise to bistability even without cooperativity between the transcription 
factors. To model this system recall that [ta] {[i^b]) can be defined as the fraction of time 
in which the promoter is occupied by a bound A (B) protein. The fraction of time in which 
the promoter is vacant is 1 — [ta] — [i^b]- Incorporating this into Eq. gives rise to the 
following modification: in the terms, each appearance of [ta] or [r^] should be replaced by 
[ta] + [^b]- For symmetric parameters, the resulting equations still exhibit a single solution, 
in which [A] = [B] and [r^] = [r^]. The Michaelis-Menten equations for the exclusive 
switch are given by Eqs. (j2)), where in the first equation k is replaced hj k/{l + k[A]) and 
in the second equation it is replaced by A;/(l + A;[-B]). To account for the discreteness of 
the transcription factors and their fluctuations, the master equation should be applied, with 
the constraint that P{Na, Nb,1,1) = 0. It takes the form of Eq. Q, except that in the 
ao and ai terms, each time 6r^j {Srsj) appears it should be multiplied by Srgfi, {^rA,o)- 
the exclusive switch, under conditions of weak repression, P{Na, Nb) exhibits a single peak 
[Fig. El^a)], for which (Na) and (Nb) coincide with [A] and [B], respectively. For strong 
repression, the distribution P{Na, Nb) exhibits two peaks. In one peak the A population 
is suppressed, while in the other peak the B population is suppressed, as expected for a 
bistable system [Fig. |2fb)]. The dead-lock situation is impossible in this system. 

To examine the time dependence of the populations of free proteins in a single cell, we 
have performed Monte Carlo simulations, based on the master equation for the exclusive 
switch. In Fig. iniwe present the copy numbers of free and bound A and B proteins vs. time. 
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The population size of the dominant specie is in the range of 20-60, while the minority specie 
is almost completely suppressed. The typical switching time is around 10^ seconds. 

Consider a different variant of the genetic switch, which exhibits bound-repressor degra- 
dation (BRD). Even a low degradation rate, rf.,., of the bound repressors tends to remove the 
mutual suppression of both species, and gives rise to a binary switch. The rate equations 
that describe this circuit are identical to Eq. (^, except that a degradation term of the 
form — ciri'^yi] {—dr[rB]) is added to the equation for [r^] ([rs]). For symmetric parameters, 
the Michaelis-Menten form of these equations, applicable in the limit of fast switching, is 
given by Eq. (j21) where k = ao/(cn + dr) and d is replaced by an effective degradation rate 
des = d + dj.k/ (1 + k[A]) in the first equation and by the analogous term in the second equa- 
tion. This equation exhibits a bifurcation at k^ = {d/dr){y/g + v^)/(v^~ in which 
the symmetric solution [A] = [B] becomes unstable, giving rise to two stable solutions in 
which one specie is dominant and the other is suppressed (Fig. HI inset). We thus find that 
in case that bound repressors exhibit degradation, bistability appears even at the level of 
rate equations. The emergence of bistability can be attributed to the fact that the effective 
degradation rate for the minority specie is larger than for the dominant specie, enhancing 
the difference between the population sizes. The master equation for this circuit is obtained 

by adding the term rf,.[(5r-,4,o - Sr,AA)Pi^A, NB,l,rB) + (5rs,o - SrB,i)PiNA, NB,rA,l)] to 
Eq. Q. This term represents transitions of the cell from rx = 1 {X = A,B) to rx = 0, 
without changing the number of free proteins. The degradation of bound repressors gives 
rise to suppression of the peak near the origin, leading to the emergence of bistability. 

A third variant of the genetic switch exhibits protein-protein interactions (PPI) such 
that an A protein and a B protein may form an AB complex, which is not active as a 
transcription factor. This circuit exhibits bistability within a range of parameters, both for 
the rate equations and for the master equation. 

We have calculated the switching time using the master equation, for an initial state that 
includes only free A proteins. The distribution P{Na, Nb) vs. time was calculated and the 
function f(t) = P{Na > Nb) — P{Na < Nb) was found to decay exponentially according 
to f{t) = exp(— t/r), where r is defined as the switching time. In Fig. 0]we present the 
switching time r, obtained from the master equation vs. k for the exclusive switch (o) and 
for the BRD switch ( x ) . We also examined the dependence of r on the copy number, A^, of 
the dominant specie. For the exclusive switch, we found that when d is varied, r ~ A^^, while 
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in case that g is varied, t ^ N. This dependence is weaker than found for the cooperative 
switch 0- 

The results presented in this paper (except for Fig. |2)) were obtained by direct integration 
of the master equation rather than by Monte Carlo methods [2^. Direct integration is 
much more efficient and provides more accurate results, without the need to accumulate 
statistics. Recent improvements in the methodology enable to use direct integration for 
complex networks that involve large numbers of active proteins , which will enable to go 
beyond elementary circuits into simulations of complete networks. 

In contrast to previous knowledge that bistability requires cooperative binding of tran- 
scription factors, we have shown that bistability is possible without cooperative binding. We 
have analyzed three variants of the genetic toggle switch, that exhibit bistability without 
cooperative binding. The first circuit is the exclusive switch, in which the two promoter 
sites cannot be occupied simultaneously. The second circuit exhibits degradation of bound 
repressors, while in the third circuit free A and B proteins may form a complex which is 
not active as a transcription factor. Rate equations predict a single stable state in the first 
circuit and bistability in the second and third circuits. However, the master equation pre- 
dicts bistability in all the three circuits. These findings are not limited to cases in which 
proteins exist in low numbers, but are due to the low copy number of the promoter itself. 
The results presented here are expected to have significant implications on the understand- 
ing of non-genetic variability in cell populations, and may shed new light on the way cells 
differentiate despite uniform environmental conditions. 

N.Q.B. was supported by the Center of Complexity of the Horowitz Foundation and the 
Bikura Program of the Israel Science Foundation. 
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FIG. 1: The probabilities P{Na, Nb) for the general switch, under conditions of (a) weak repression 
(fc = 0.005) where there is one symmetric peak and (b) strong repression (k = 50) where three 
peaks appear, one dominated by A, the second dominated by B and the third in which both species 
are mutually suppressed. The weights of the three peaks are about the same. 

FIG. 2: The probabilities P{Na,Nb) for the exclusive switch, under conditions of (a) weak 
repression {k = 0.005) where there is one symmetric peak and (b) strong repression {k = 50) where 
bistability is observed. 

FIG. 3: The populations of free and bound A and B proteins vs. time, obtained from Monte 
Carlo simulations of the exclusive switch with the parameters g = 0.2, d = 0.005, ag = 0.2 and 
ai = 0.01. The bistable behavior is clearly observed, where the population size of the dominant 
specie is between 20-60 and the other specie is nearly diminished. Failed switching attempts are 
clearly seen. 

FIG. 4: The switching time vs. the repression strength, k, for the exclusive switch (o) and for the 
case in which bound repressors exhibit degradation (x). For the bistable range (roughly k > 1) the 
switching time increases as k is increased. The inset shows the steady state solution for [A] and 
[B] vs. k, obtained from the rate equations for the BRD switch. Note that for the BRD switch, 
the parameter oq varies, while dr = d and ai = 0.01 are held fixed. 
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